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Abstract 

It is shown that the one- and two-dimensional solids may exist and 
melt. 

Suppose that we have N identical atoms, of mass M each, arranged 
along a rectilinear chain; suppose that the atoms interact through a nearest- 
neighbours potential, and we are interested in small atomic displacements 
Ui , i = 1, 2, ...N, of immaterial polarization, around equidistant equilibrium 
positions Xi — i ■ a, where a is the lattice constant. Within the harmonic 
approximation with the elastic force constant K this motion is described by 
the hamiltonian 

# = E^ 2 + ^£K-%) 2 . a) 

i (ij) 

where pi is the conjugate momentum, [pi, Uj] = —ihSij . With cyclic boundary 
conditions we have 



Ui = 




(2) 



k = (2n/N a) -integer, —n/a < k < ir/a, u\ = , 



Pi = M Ui= ~y= ^Zpke ll 




(3) 
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p k = M u_ k , [p k , u k >] = —ih5 kk ' , and the hamiltonian (1) becomes 

H = E {^PkPk + lMu 2 k u{ Uk ] , (4) 

k 

where 

AK 

w * = w^ 1 ~ coska "> ; ( 5 ) 

it is brought into the diagonal form 

# = 5>u; fc (n fe + 1/2) , (6) 
k 

where n k = ala k , by introducing the phonon operators a k , 



^JL-(a k + al h ) , Ph = i^™pL(al-a- k ) . (7) 

The mean square deviation of the atomic displacements in the phonon vac- 
uum 

diverges logarithmically at k — . This infrared divergence, in various con- 
texts, is currently invoked to rule out the existence of the one- dimensional 
solid. i'§ 

However, the mean square deviation of the lattice constant 

L„\ — 

(9) 



= («j - Mj) = jjY,k u k u k -2(1- cos ka) = 
& £ fc jjfer (1 - cos feo) (2n, + 1) = f J^ dk-jt-Al- cos fca) 



is finite, which indicates that the solid exists, only its center of mass fluctuates 
immeasurably. Obviously, the thermodynamics is meaningless in this case. 

The question one should ask is whether the phonons do exist in the one- 
dimensional solid subjected to the actual experimental conditions. These 
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always imply an uncertainty u in the localization of the center of mass of 
the rigid lattice, and the constraint under which we should look for phonons 
reads 

J2uf = Nu* . (10) 

i 

The magnitude of u can not be smaller than a/2 , unless the crystal is 
destroyed, but usually it is much larger than the lattice constant a. 



Introducing the Lagrange multiplier oj^ and setting 

H = H-\Mu J lY J {u 2 -ui) (11) 
we get straightworwardly 

w = E (2^ p£p* + ^ Mn Wk) - ^o 2 ^ 2 , (12) 

where 

^ = ^ 2 + ^ , (13) 

z. e. the phonon spectrum acquires a threshold frequency u which acts as an 
infrared cut-off. In addition, the ground-state energy 

E = H = J2 \*tok - \mujINu 2 (14) 

k 

has a minimum with respect to Uq, 



dE _ „ ft 1 B/f . Ar __ 2 
~ V 4&h ~ 2 



which is exactly the constraint (10) ; indeed, 



E^->^ 2 , as) 



ft 



ik k h 

coincides with (15). Assuming a Debye spectrum u K = ck , where c = 
a(2K I 'M)a is the sound velocity, and introducing the Debye wavevector fc^i = 
7r/a and the Debye frequency ujd = ckn , we get from (16) 

VD/U ° dx- 1 - (17) 



/ 

Jo 



The right-hand side of (17) is, typically, very large, so that we get the thresh- 
old frequency 

2Mwn n 

u = 2u D -e —' u (18) 

which is extremely small. 



The phonon partition function for 7i given by (12) reads 

Z = TT — = e -( F +k M "l N ^)/T fjn) 

^ lll_ e -nn k /T e ' ^ iy J 



k 



where F is the free energy; its minimum dF/ du^ = gives 

i.e. again the constraint (10) , at finite temperatures. Equation (20) reads 

r°-7r^- h (^)^'" 2 ' (21) 

whose solution for T — > is (18) and 

for T — > oo ; this is again an extremely small frequency for any normal 
temperature, and (22) may be used satisfactorily for any finite temperature. 
The constraint (10) introduces therefore an extremely small infrared cut- 
off ujq in the phonon spectrum, according to (13), slightly depending on 
temperature, whose effect is practically unobservable. We remark that the 
mean square deviation of the distance between any pair of atoms 

h 



(ui - Uj f = -J2 jTTTj- ' 2 {1 - cos [ka (< - j)]} ■ (2n k + 1) (23) 

k k 

is always smaller than 4-u 2 , according to (20) , so that the (on-diagonal) 
crystalline long-range order does exist in one dimension. 
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In this connection the question of melting of the one- dimensional solid is 
meaningful! The elastic force constant K in (1) is a function of a ; actually 
it is a function of the nearest-neighbours distance a + Ui — Uj ; more than 
this, it is a periodic function with the period a , whose integral extended to 
the whole crystal vanishes. Consequently, it may be expanded as a Fourier 
series of the reciprocal vectors of the lattice, and we may write it as 

K(a + Ui — Uj) = K cos [G (uj — Uj)\ + ... , (24) 

where G = 2ir/a . Assuming vanishing fluctuations we may approximate the 
cosine in (24) by its average 



cos[G(u, -Uj)} = e"^ 2 , (25) 



where 5 2 = [G (ui — Uj)] 2 , and keep the first G-term in (24). The sound 
velocity changes, therefore, from c to c • exp (— 5 2 /4) , and we get a mean- 
field theory known as the self-consistent harmonic approximation. Remark 
that it is (ui — Uj) 2 which enters the theory, not uj . On the other hand 

5 2 = G 2 (u t - Uj) 2 = jtY. ( X ~ cos ka ) ■ (2** + !) > ( 26 ) 

k k 

and we may neglect Uq here and approximate 1 — cos ka by (ka) 2 /2 over the 
whole integration range. With the Debye spectrum (26) becomes 

5 2 = • T 2 ■ et- 52 ■ " dx • xcothx . (27) 

hMu% Jo v 1 

Threre is a critical temperature beyond which this equation has no solution 
anymore; it is given by the tangent point of the curves described by the two 
sides of (27) , which may be approximated by 



hMu 3 n ' ' 4 



whence, together with (27), yields 



D, 



■ T 2 ■ -e* 5 ■ / dx-xcothx , (28) 



«=j- (29) 
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The high-temperature limit of (27) reads now 



7T 2 G 2 



f = ^—. T -e^ 2 , (30) 



whence we get the melting temperature (not to be mistaken for the freezing 
temperature §]) 



It is interesting to compare (u{ — Uj) 2 = a 2 /3ir 2 at the melting point given 
by (29) with the square of the lattice constant a 2 . The melting temperature 
given by (31) indicates a sharp transition where the crystal gets soft and 
can no longer bear phonons. The investigation of the validity of this mean- 
field theory may indicate, in general, a continuous transition, with a variable 
range of crystallinity, both in magnitude and orientation. H Indeed, the well- 
known [|] argument for the existence of the multi-phased one- dimensional 
thermodynamical systems can be applied here for the breaking up of the inter- 
atomic bonds between neighbouring atoms. For a one-dimensional crystalline 
solid consisting of N atoms, with an energy if) for each bond, the equilibrium 
is reached for a certain, finite value of the concentration n/iV C 1 given by 



n 



N-e~^ T , (32) 



where n is the total number of broken bonds. This shows that the one- 
dimensional crystalline solid is unstable, at any finite temperature, with re- 
spect of the breaking up of the interatomic bonds (including the bonds with 
the substrate that ensure the constraint given by (10)). This instability may 
be viewed as a continuous transition toward a state with a variable range of 
crystallinity. We should remark, however, that the concentration given by 
(32) is extremely small for any realistic values of the bonding energy if) and 
normal temperatures. 



The constraint (10) applies to all the atoms in the lattice, as, for example, 
to a lattice with a basis in the unit cell. However, the optical phonons do not 
contribute practically to the translational displacement of the lattice, so that 
the threshold frequency ujq (appearing now in the whole phonon spectrum) 
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is practically unchanged. The constraint (10) applies also to solids of any 
dimensionality, and we shall discuss now, for the sake of reference, the case 
of a cubic three-dimensional Bravais lattice. The hamiltonian of the atomic 
vibrations can be written as 

H = Y.J^Pli + \ 2 G «P ( r * - r i) u ™ u Pi > ( 33 ) 

cti oc/3{ij) 

where a, f3 are the cartesian labels of the components of the displacement 
vector Uj , are the equilibrium positions, and G a p is the matrix of the 
elastic force constants. A Fourier expansion brings H into 



1 1 

H = H —A ■ PlkPak + 77 S G ^ W ' M «k M /3k (34) 

where 



2M " ak "" K ' 2 „ 

ak apk 



G aP (k) = ]T Gap (r, - r,-) • , (35) 

j 

j's being the nearest neighbours of % . A canonical transform diagonalizes 
G a/ 3 (k) into G a (k) , where we use the same label a for polarizations. In 
addition, we assume a Debye, isotropic, identical spectrum for all the polar- 
izations, 

G a (k) = Mul = Mc 2 k 2 . (36) 
We remark that the constraint (10), which now reads 

£^ = A M 2 (37) 

i 

is left unchanged under the canonical transform. Under these conditions we 
get the spectrum VL 2 k = oj 2 + uj\ and (37) gives 

E ~4tt ■ ™th ( = Nu 2 . (38) 

1 /3 

Making use of the Debye wavevector hp = (Qir 2 ) /a , and the Debye fre- 
quency ujd = eke , (38) becomes 

SHU 2 fUD/uJo x 2 / % Uo \ u 2 



dx-^=L=.coth r^y/T+&\=±. . (39) 



2Mu% Jo VT+x 2 V 2T 



For T — > the only solution of (39) is uoq = and the zero-point fluctuations 
Uq = 9h/4MujD (<C a 2 ) ; for T — > oo we find wo = again, and w 2 = 
9T /Muj 2 d , which varies extremely slowly with the temperature, and may be 
used for any normal range of temperatures. Therefore, there is no frequency 
cut-off on the phonon spectrum in three dimensions, and the lattice develops 
its own fluctuations. 



In order to compute the melting temperature we use again the self- 
consistent harmonic approximation, according to which 



G a(3 (Ti -Tj)=> G a(3 (Ti - Tj + Ui- Uj 



Eg Gap (G) cos [G (u, - u,)] Eg G a p (G) e -|[G(ui-»i)] a 
On the other hand 

G (llj - Uj) = E« G a (U ai - U a j) = 



(40) 



and 



G fa - uj)] 2 = ^ T Y. G l- M lk M «k • 2 {1 - cos [k (n - rj)]} 



,ikri 



N 



(41) 



(42) 



ak 



For a cubic lattice the sum in (42) does not depend on (ij) and it may be 
approximated by 



[G(u,-u 3 )r 



G 2 a 2 
N 



E 



2Mu k 



k 2 x -(2n k + l) , 



(43) 



or 



[G (Ui - u,)] 2 = G 2 8 ^ 1 ' • T 4 r DlZ± dx ■ x 3 cothrr . (44) 
n rMujp Jo 

We may therefore keep the first G s in (40) and, denoting the corresponding 
(44) by 5 2 , we have 



2\2/3 



hu D /2T 



G^ {ji - Tj) =>• G a/3 {Ti - Yj) ■ e 2 



5 2 



-±<5 2 

c • e 4° 



(45) 
(46) 
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and a similar relation for ujd ■ Equation (44) becomes 

ft Mwf, jo 
whose critical solutions are given by 

<5 2 =4/5 . (48) 
The finite-temperature estimation of (47) is 

<* = IS , («) 

whence the melting temperature 

T m = ■ M ° 2 - 0M ■ M ° 2 ; ( 5 °) 

this is the Lindemann's law of melting. || At the melting temperature the 
fluctuations in the atomic positions are 



9 

>D 



U ™ = ' Tm ~ °'° 25a2 ' (51) 



according to (39) , and we may see that their ratio to the zero-point fluctu- 
ations is given by 

u m/ u o = 4T m /huj£) = 0.04 • ■ (52) 



We shall discuss now briefly the two-dimensional case. Under the same 
assumptions as those employed for the three-dimensional solid the constraint 
equation (39) becomes 

^ ^'"V.^.cothf^v/TT^U^ , (53) 



Mul Jo VTTx* V 2T 

where ujd = ck^, = 2^fn/a . For T — > this equation has the solution 
ujq = and the zero-point fluctuations Uq = 2h/Muj£,. For T — > oo we get 

u = u D - e~—- u , (54) 
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i.e. the two-dimensional solid may exist at finite temperatures, [0 under 
the constraint (10). The cut-off frequency given by (54) is extremely small 
and varies very slowly with the temperature; it may be used for any finite 
temperature range (remark that it vanishes for T — » ), and its effects 
are practically unnoticeable. The melting may be treated similarly with the 
three-dimensional case. Instead of (47) we get now 

52 = G 2 -2^— - • T 3 • e 52 *' dx ■ x 2 ■ cothx , (55) 

n Muf) Jo 

whose critical solution is given by 

5 2 ^ 1 . (56) 
The finite-temperature integral in (55) gives 

2 = 2^ 

Mc 2 

whence the melting temperature 



T-e* s , (57) 



T ™ = ^tte ■ M ° 2 ~ °- 03 • M ° 2 ■ ( 58 ) 

We stress again the possibility of a continuous transition, |Q including various 
structural phases. 



Similar considerations as those presented here apply for highly anisotropic 
solids. A natural question arises of what happens for quasi-low dimensional 
solids, as for a slab or a rod of thickness d . The qualitative picture given here 
remains unchanged, but, of course, the quantitative results are altered. The 
main feature is the occurrence of a cross-over temperature toward a three- 
dimensional behaviour, of the order of hc/d ; for d = N t a this temperature 
is about 10 2 /N t Kelvins, for a typical sound velocity c ~ 10 3 m/s . 
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